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ABSTRACT 


When  the  predicted  position  of  a  satellite  contains  normally 
distributed  errors,  the  position  uncertainty  can  be  described  by  a 
Spherical  Error  Probable  or  SEP.  The  SEP  is  calculated  by  inte¬ 
grating  the  three-dimensional  normal  probability  density  function 
over  a  spherical  volume.  The  SEP  is  set  equal  to  the  radius  of 
that  volume  which  contains  the  satellite  with  50%  probability.  In 
this  report  the  authors  present  four  methods  for  integrating  the 
density  function  and  finding  the  SEP.  The  three  normal  variates 
in  the  density  function  are  assumed  to  be  independent  and  unbiased 


with  known  variances. 
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SECTION  1 


INTRODUCTION 


1.  1  Purpose  of  the  Report 

In  this  paper  we  report  four  methods  of  solving  Eq.  (1-1)  for  R 
when  ando_  are  known  and  P  is  set  equal  to  0.  5. 

Jv  V  2 


dx  dy  dz 


d-1) 


Eq.  (1-1)  describes  the  integration  of  an  unbiased,  three-dimensional, 
normal  probability  density  function  over  a  spherical  volume.  R  is  the 
radius  of  the  sphere;  x,  y  and  z  are  independent  random  variables, 
and  P  is  a  probability.  For  P  set  equal  to  0.  5,  we  define  R  as  the 
Spherical  Error  Probable  (SEP). 

I 

1.  2  Uses  for  the  SEP 

The  SEP  is  useful  as  a  measure  of  position  uncertainty  in  three 
dimensions.  It  is  a  single  number,  like  3250  ft,  describing  the  un¬ 
certainty  in  the  location  of  one  vehicle  at  a  given  time.  The  SEP  is 
applicable  to  problems  in  satellite  tracking,  missile  defense,  sub¬ 
marine  navigation  and  air  traffic  control.  ^)^)(7)  As  an  example, 


I 


we  describe  in  the  next  four  paragraphs  a  SEP  that  might  be  computed 

in  predicting  the  location  of  a  manned  orbiting  laboratory. 

.  • 

A  manned  orbiting  laboratory  in  a  near-earth  orbit  requires  re- 
supply.  The  supply  mission  is  to  be  flown  by  a  second  missile,  the 
transporter,  and  a  rendezvous  and  docking  maneuver  is  to  be  executed 
when  the  transporter  reaches  the  laboratory.  A  successful  mission 
hinges  on  an  accurate  prediction  of  the  laboratory's  location  at  the 
time  of  rendezvous.  The  mission  planners  are  asked  to  predict  the 
laboratory's  future  position  by  using  radar  measurements  made  just 
before  the  transporter  is  launched. 

The  mission  planners  would  like  to  make  an  absolute  prediction 
like  the  following: 

"The  laboratory  will  be  at  latitude  20°  N,  longitude  35°  W 
and  100  nm  altitude  at  2117  Greenwich  Mean  Time.  " 

They  cannot  make  such  a  statement,  however,  because  the  radar  data 
contains  normally  distributed  errors.  The  present  location  of  the  lab- 
oratqry  is  uncertain  so  its  future  position  cannot  be  predicted  exactly; 
the  predicted  latitude,  longitude  and  altitude  contain  the  unknown  ran¬ 
dom  errors  x,  y  and  z,  respectively. 

The  mission  planners  choose  to  compute  a  SEP  to  describe  the  un¬ 
certainty  in  the  future  position  of  the  laboratory.  Their  computed  SEP 
iis  8000  feet.  The  mission  planners  use  the  SEP  to  make  the  following 


probablistic  statement  to  the  transporter  crew: 

"The  laboratory  will  be  within  8000  feet  of  20°  N 

latitude,  35°  W  longitude  and  100  nm  altitude  at  2117 

Greenwich  Mean  Time  with  50%  probability." 

The  transporter  crew  thus  has  an  estimate  for  the  uncertainty  in  the 
laboratory's  location.  Upon  reaching  the  predicted  point  of  rendez¬ 
vous,  the  crew  will  probably  have  to  search  through  a  sphere  of  at 
least  8000  foot  radius  to  find  the  laboratory. 

The  previous  example  illustrates  the  information  implied  in  the 
Spherical  Error  Probable.  If  a  vehicle  is  said  to  be  located  at  point 
L,  K  and  H  with  a  SEP  of  8000  feet,  we  immediately  know  that: 

(1)  The  location  of  the  vehicle  is  uncertain. 

(2)  The  three  position  coordinates,  L,  K  and  H,  contain  normally 
distributed  errors. 

(3)  The  probability  density  function  is  so  shaped  that  50%  of  the 
probability  is  contained  in  a  sphere  of  8000  foot  radius.  The 
sphere  is  centered  on  the  mean  predicted  location  of  the 
vehicle. 

A  SEP  can  be  formulated  to  describe  position  uncertainties  for 
many  vehicles  -  satellites,  aircraft,  submarines,  and  missiles.  In 
this  report,  however,  we  restrict  the  SEP  application  to  those  three- 
dimensional  problems  where  the  random  position  errors  are  normally 
distributed,  unbiased  and  independent. 
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1. 3  The  ‘Three-Dimensional  Normal  Distribution 
- - 1 -  - 

If  three  random  variables,  x,  y  and  z,  are  independent,  unbiased 
and  normally  distributed,  their  joint  three-dimensional  probability 
density  function  is  given  by  the  equation^: 


f  (x,  y,  z) 


1 _ 

(2^ /2o  o  r 

x  y  z 


exp 


If  the  density  function,  f  (x,  y,  z),  is  integrated  over  a  closed  volume, 
a  probability,  P,  results.  If  the  three  random  variables  describe  the 
position  of  a  point,  P  is  the  probability  that  the  point  is  located  some¬ 
where  in  the  volume.  In  this  report  the  closed  volume  is  a  sphere. 

R  is  defined  as  the  radius  of  the  sphere  and  Eq.  (1-1)  describes  the 

(2) 

probability  integral  for  Px  ;  ^ 

r  Vr2-^-/ 


When  the  standard  deviations,  o  ,  a  o  ,  are  known  and  R  is 

^  /  Z 

( l )  (2 ) 

specified,  Eq.  (1-1)  can  be  solved  for  P.  We,  however,  are 

interested  in  the  complementary  problem.  We  seek  a  solution  for  R 


when  P  is  given.  In  particular  we  want  R  wheno  o  and  o  are 

*•  V  Z 


4 


known  and  P  is  set  equal  to  0.  5.  This  R  is  defined  as  the  Spherical 
Error  Probable  or  the  SEP. 

The  probability  limit,  in  this  case  0.  5  or  50%,  is  not  unique. 

Any  limit,  say  95%,  can  be  used  in  the  definition  of  the  SEP.  The 
50%  limit,  however,  is  consistent  with  the  definition  of  the  two-dimen¬ 
sional  Circular  Error  Probable  (CEP)  used  in  missile  accuracy 
studies.  ^ 

1 . 4  Outline  of  the  Paper 

In  the  next  four  sections  we  outline  four  methods  for  computing 
the  SEP  when0  ,  a  and  o  are  given.  In  Section  2,  for  example,  we 

X  y  Z 

use  the  CEP  curve  to  approximate  the  SEP  in  the  special  case  when 
0  =  o  in  Section  3  a  paper  and  pencil  solution  is  developed. 

x  y 

Section  4  describes  a  computer  method  for  finding  the  SEP,  and 
Section  5  summarizes  the  computer  data  in  a  graphical  solution. 

The  paper  closes  with  three  appendices.  The  first  appendix 
extends  the  approximation  method  of  Section  2  to  the  cases  where 

I 

°x  4  °y  The  second  appendix  outlines  H.  W.  Lilliefors'  solution 
of  Eq.  (1-1).  The  third  appendix  contains  a  computer  program  and 
sample  data. 

The  reader  with  standard  deviations,  °  ,  0  and  0  ,  in  hand  and 

x  y  z 

a  deadline  to  meet  should  proceed  directly  to  Fig.  5.2  in  Section  5. 
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For  a  wide  range  of  standard  deviations,  Fig.  5.2  permits  a  direct 
reading  of  the  SEP. 
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SECTION  2 


THE  FIRST  SOLUTION:  A  METHOD  OF  ESTIMATING  THE  SEP 

2.  1  Outline  of  the  Estimation  Method 

In  general,  the  SEP  cannot  be  found  by  a  direct  integration  of  Eq. 

(1-1).  We  cannot  write  a  mathematical  equation  that  expresses  R  in 

terms  of  P,  o  ,  o  and  o  .  We  can,  however,  estimate  the  SEP  by 
x  y  z 

looking  at  special  cases.  In  this  section  we  show  how  the  SEP  can  be 
estimated  when  o  and  a  are  equal  and 

X  y 

(1)  Small  compared  to  oz. 

(2)  Large  compared  to  a  z. 

(3)  Equal  to  oz. 

(4)  Equal  to  one -half  of  o 

Here  we  also  introduce  the  normalized  function,  SEP/oz,  because  it 

is  more  easily  plotted  than  the  SEP. 

2.2  The  SEP  When  o  and  o  are  Equal  and  Small 
i _ x  °y _ 2 _ 

In  Sec.  1.  3  we  wrote  down  Eq.  (1-2),  the  three-dimensional 
normal  density  function: 


f  (x,  y,  z)  = 


1 


/?_v3/2 

(2’)  °X°yOz 


exp 


-*(£*:£♦£)]  <« 


If  two  of  the  standard  deviations  in  this  expression  are  equal  and 


smajiler  than  the  third  deviation,  the  three-dimensional  distribution 


look#  one-dimensional.  For  example,  if  ox  and  ay  are  equal  and 
smaller  than  0Z,  the  three-dimensional  distribution  looks  like  a  one* 
dimensional  distribution  in  the  variate  z.  In  the  limit,  as  the  ratio 
approaches  zero,  the  SEP  calculation  reduces  to  the  problem 
of  finding  the  half-length  of  a  straight  line.  The  full-length  of  the 
straight  line  includes  the  location  of  the  random  point  with  50%  prob 
ability. 

The  probability  equation  to  be  solved  is  just: 


(2-1) 


We  seek  R  when  o  is  given  and  P  is  set  equal  to  0.  5. 

Z 

If  q  is  substituted  for  z/oz,  the  integrand  in  Eq.  (2-1)  becomes 

the  standard  normal  density  function  for  which  there  are  tabulated 

(4) 

solutions.  Writing  Eq.  (2-1)  in  the  standard  form  we  get: 


(2-2) 


-R/°» 

We  use  the  tables  and  find  that  when  R/oz  =  0.674,  P  equals  0.  5. 
Thus,  for<*x  and  <*y  equal  to  zero,  the  ratio  SEP/°Z  has  the  value 
0.  674.  This  result  is  plotted  in  Figure  2.  1,  Sec.  2.  6,  as  a  point 


on  the  ordinate  where  ox/az  =  0.  This  point  is  plotted  as  the  first 
step  in  constructing  a  curve  of  SEP/^Z  versus  ox/o z  under  the  special 
condition  that  ox  equals  a  y. 

2.  3  The  SEP  Whenox  and  uv  are  Equal  and  Large 

We  have  seen  that  the  three-dimensional  density  function  looks 
one-dimensional  when  two  of  the  standard  deviations  are  equal  and 
small  compared  to  the  third.  If,  conversely,  two  of  the  standard 
deviations  are  equal  and  large  compared  to  the  third  deviation,  the 
distribution  looks  two-dimensional.  If  ox  and  o  are  equal  and  much 
larger  than  az,  the  dens  it/  function  looks  two-dimensional  in  the 
variates  x  and  y. 

Calculation  of  the  SEP  reduces  to  the  problem  of  finding  the  radius 
of  a  circle.  The  circular  area  contains  the  location  of  the  random 
point  with  50%  probability.  The  desired  radius  is  commonly  called 
the  Circular  Error  Probable  (CEP).  The  CEP  is  usually  used  to 
specify  the  probable  impact  error  of  a  long  range  missile. 

The  two-dimensional  probability  integral  is  given  by  Eq.  (2*3). 


R  tJfP-J- 
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A  solution  for  this  integral  is  graphed  in  Ref.  3.  For  0y/ax  greater 

i 

than  1/3  the  SEP  is  approximated  by  Eq.  (2- 4). 


SEP  =  0.  59  (°x  +  o  ) 


(2-4) 


In  our  case,  with  Oy  and  ox  equal  and  much  larger  thanc^,  the  SEP  is 


given  by  the  expression 


SEP  =  1.  18  °x  =  1.  18  oy 


(2-5) 


Equation  (2-5)  is  made  more  convenient  if  we  divide  both  sides  of 
the  equation  by  o  . 


^  =  1.18  ^ 


(2-6) 


This  result  is  plotted  in  Fig.  2.  1  as  a  straight  line  with  slope  equal 
to  1.  18  and  intercept  at  c»x/  oz  =  0.  The  curve  for  ox  =  oy ,  in  Fig. 
5.2,  Sec.  5,  can  be  used  to  verify  that  this  expression  for  the  SEP 
is  accurate  when  the  ratio  ax/oz  is  greater  than  approximately  4.  5. 
2.4  The  SEP  when  ox,  Qy  and  are  All  Equal 

iWhen  all  three  standard  deviations,  ox,  a y  andoz,  are  equal,  the 

i 

probability  density  function  retains  its  three-dimensional  form.  We, 

however,  can  simplify  the  integral  equation  for  P  by  changing  to 

2  2  2  2 

spherical  coordinates  and  substituting  P  for  x  +  y  U  in  the  inte¬ 


grand.  The  new  expression  is  then: 


10 


(2Ti  o. 


R  2*  it/  2 

Jin 

0  0  -n/2 


j  P  sin$dpd0d*  (2-7) 


J.  S.  Toma,  in  Ref.  1,  reports  that  if  P  =  0.  5,  Eq.  (2-7)  can  be 


solved  to  give  the  result  that  the 


SEP  =  1.53820 


(2-8) 


If  we  divide  both  sides  of  this  expression  by  0  ,  Eq.  (2-9)  appears: 

z 


SEP 


=  1.5382 


(2-9) 


Thus,  for  °x  =  °y  -  az  we  have  the  result  that  SEP/  °z  =  !•  5382.  This 
result  is  plotted  in  Fig.  2.  1  as  a  point  at  °x/°z  =  1. 

2.5  The  SEP  when  qx  =  °y  =  0.  5  <*z 

The  last  special  case  we  examine  is  that  for  o  and  o  equal  to 

x  y 

0.  5 o  The  SEP  is  calculated  using  the  work  of  H.  VV.  Lilliefors  in 
z 

Ref.  2.  He  plots  the  probability,  P,  versus  o ^  when  R  =  1.  His  Fig. 
2  shows  that  P  :  0.5  when  °x  =  0.  5,  °y  =  0.  5  and  °z  =  1.  Thus,  for 
ox/oz  =  0.  5,  the  ratio,  SEP/oz>  must  be  equal  to  one. 

2.6  The  Predicted  SEP  Curve  for  o  and  °v  Equal 


In  Sec.  2.2  through  2.  5  we  have  computed  one  line  and  three  dis¬ 
crete  points  which  describe  the  behavior  of  the  function  SEP/°Z.  We 
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have  plotted  these  points  and  the  line  in  Fig.  2.  1  and  we  now  draw  in 


a  dotted  line  to  complete  the  curve.  This  graph  can  be  used  to  esti- 


♦ 


1  2  3  4  5  6  7 


a 

x 


Figure  2.  1  The  predicted  SEP/&Z  curve  when  o  x  =  oy. 


special  condition.  When  ox  and  are  not  equal,  the  SEP  must  be 
estimated  by  the  methods  of  Appendix  A. 


*  The  curve  can  be  used  to  comp  '.te  the  SEP  when  any  two  stand¬ 
ard  deviations  are  equal.  The  reader  need  only  interchange  the 
subscripts,  x,  y  and  z,  on  the  standard  deviations. 


To  find  a  SEP  given  ox  =  oy,  oz  and  Fig.  2.  1,  we  first  find  the 
ratio  °x/az.  If  the  ratio  has  the  value  2.2  5,  for  example,  we  enter 
Fig.  2.  1  at  the  bottom  where  ox/oz  =  2.25.  We  proceed  vertically 
upward  to  the  curve  and  read  that  the  SEP/oz  has  the  value  2.85. 

The  SEP  is  then  computed  by  multiplying  2  85  by  the  standard  devia¬ 
tion,  °z. 

A  more  exact  solution  can  be  found  by  using  paper,  pencil,  a  desk 
calculator  and  the  iterative  method  of  Sec.  3. 
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SECTION  3 


THE  SECOND  SOLUTION: 

AN  ITERATIVE  METHOD  OF  FINDING  THE  SEP 


3.  1  Outline  of  the  Method 

Given  one  set  of  standard  deviations  and  P  set  equal  to  0.  5,  a 
SEP  can  be  calculated  using  paper,  pencil  and  a  desk  calculator.  The 
method  uses  H.  W.  Lilliefors'  solution  for  Eq.  (1-1)  and  we  solve  by 
iteration  for  the  SEP. 

3.  2  Lilliefors1  Solution 

Lilliefors'  solution  (see  Appendix  B)  for  the  three  dimensional 


probability,  P,  has  the  form: 


(1) 


.ysT?  f  (-i)(ntl)  tt|<c^ n)  - 


1 


(j+n+j)  (j+n-j).  •  •  (j+j) 


(-i  r 


j!(2az2)j 


(3-1) 


This  solution  is  valid  when  R  =  1.  For  R  not  equal  to  one,  Lilliefors1 


solution  has  the  revised  form: 


p  =  yfzh  £ 


(n+l) 

(-1)  «!(C2,n) 


n«l  2 


(n+l)  .  v(2n-l) 


C4  f*)f+) 


j-0 


(-1) 


0 +n+})  (j+n-I)...(j+i) 


(3-2) 
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To  find  the  SEP,  we  set  P  equal  to  0.  5  and  solve  for  R.  The  equation, 
however,  is  too  complicated  to  be  solved  in  closed  form.  We  use  an 
iterative  procedure. 

3.  3  The  Iterative  Procedure 

The  first  step  in  the  iteration  is  to  specify  °x,  °y  and  °z<  Next, 
we  make  an  initial  estimate  fore  the  SEP,  say  R  =  RQ.  ^0,  ox,  Oy  and 
°z  are  then  substituted  into  Eq.  (3-2)  and  we  solve  for  P.  If  each  of 
the  normalized  standard  deviations,  0x/Ro.  0y/RQ  an<*  °Z/R0,  is 
greater  than  about  0.4,  the  two  infinite  series  in  Eq.  (3-2)  can  each 
be  truncated  after  10  terms. 

The  solution  proceeds  by  comparing  the  computed  P,  say  it's  PQ, 
with  0.  5.  If  PQ  is  equal  to  0.  5,  we  define  the  initial  estimate,  RQ,  as 
the  SEP.  The  solution  is  deemed  complete  for  the  given  values  of  °x, 
Oy  and  oz.  If,  however,  Po  is  not  equal  to  0.  5,  we  make  a  new  esti¬ 
mate  for  the  SEP,  say  R  =  Rj,  and  repeat  the  procedure. 

If  the  reader  has  one  set  of  standard  deviations,  paper,  pencil  and 

i 

some  patience,  this  procedure  will  produce  a  SEP  in  a  few  hours.  The 
calculation  time  is  reduced  if  you  have  a  good  first  estimate  for  the 
SEP.  If,  however,  many  sets  of  standard  deviations  are  to  be  used, 
the  time  required  to  calculate  the  exact  SEPs  is  enormous.  In  this 
case  we  recommend  the  computer  solution  in  Section  4. 
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SECTION  4 


THE  THIRD  SOLUTION: 

A  COMPUTER  METHOD  OF  FINDING  THE  SEP 


4.  1  Outline  of  the  Method 

The  computer  solution  for  the  SEP  mechanizes  the  iterative  pro- 
cedure  outlined  in  Sec.  3.  The  computer  method,  however,  forces 
us  to  deal  with  three  mathematical  problems  that  we  ignored  in  Sec.  3. 
First,  we  have  to  identify  the  sets  of  standard  deviations  for  which 
Lilliefors'  series  solution  will  not  converge.  Second,  we  must  find 
alternate  methods  for  finding  the  SEP  when  the  series  solution  diverges. 
Third,  we  need  to  write  a  mathematical  rule  for  estimating  the  SEP. 

4.  2  Conditions  for  Convergence 

H.  W.  Lilliefors'  series  solution  for  P  is  used  in  the  search  for 
the  SEP.  The  revised  form  of  Lilliefors'  solution  is  given  by  Eq.  (4-1). 


(j+n+j)  (j+n- j).  .  .  (j+j) 


(4-1) 


*  A  listing  of  the  computer  solution  is  presented  in  Appendix  C. 


When  any  one  of  the  normalized  standard  deviations,  o  /R,  o  /R  or 

x  y 

°X/R,  is  less  than  approximately  0.2,  Eq.  (4-1)  will  not  converge  to 

* 

a  meaningful  value  for  P.  If,  for  example,  °Z/SEP  =  0.11  and  100 
terms  are  used  in  each  series,  Eq.  (4-1)  may  give  a  probability  like 
2  x  10^,  or  -0.93,  or  25.7.  None  of  these  numbers  are  admissible 
values  for  P.  The  probability  can  only  have  a  value  between  zero  and 
one.  When  calculating  the  SEPs  for  many  sets  of  standard  deviations, 
we  must  identify  those  sets  of  standard  deviations  for  which  the  com¬ 
puter  solution  will  not  converge. 

Inadmissible  sets  of  standard  deviations  are  identified  by  testing 

the  ratios  0  /SEP,  °V/SEP  and  0 '/SEP  against  0.2,  the  convergence 
x  y  z  “ 

limit.  We  form  the  ratios  by  dividing  the  standard  deviations  by  an 
estimate  for  the  SEP.  The  SEP  is  estimated  by  the  methods  of  Sec.  2 
or  Appendix  A.  If  one  of  the  ratios,  °Z/SEP  for  example,  is  less  than 
0.2,  we  conclude  that  Eq.  (4-1)  cannot  be  used  to  find  the  exact  SEP 
for  that  set  of  standard  deviations.  We  look  for  another  method  of 
finding  the  SEP  for  that  set. 


*  We  are  not  concerned  here  with  absolute  convergence  or 
divergence.  We  define  Lilliefors'  solution  to  be  convergent 
if  we  get  a  meaningful  value  for  P  by  using  a  reasonable 
number  of  terms  in  each  of  the  two  infinite  series.  In  this 
paper  we  call  100  terms  reasonable. 
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4.  3  Alternate  Solutions  for  the  SEP 


If  a  set  of  three  deviations  contains  one  element,  for  example, 

Z 


that  fails  the  convergence  test,  we  can  choose  to  ignore  o  and  solve 

z 


for  a  SEP  in  two  dimensions.  The  SEP,  so  calculated,  is  obviously 


an  approximation  to  the  true  SEP  since  we  assume  that  °L  =  0  when, 


in  fact,  it  is  a  non-zero  positive  number.  The  approximation  is  good. 


however,  when  0„/SEP  is  less  than  about  0.2. 


The  two-dimensional  SEP  can  be  calculated  directly  from  Eq. 

(3) 

\  i _ _ _ i:  _  n  in  •_ _ . _ .i _ _  i _ . _ j _ i  \  > 


(4-2)  when  the  ratio,  °y /  °x,  is  greater  than  about  one -third. 


SEP  =  0.  59  (°x  +  o  ) 


(4-2) 


The  SEP  can  also  be  calculated  by  iteration  if  we  substitute  Eq.  (4-3) 


for  Eq.  (4-1)  in  the  computer  solution. 


p  =  y;  — 

on 


<-l)(n+1)  «i<C2,  n) 


2"  n!  (0y/R)(2"7T>('’x/R) 


(4-3) 


The  computer  program  in  Appendix  C  incorporates  both  equations 


(4-1)  and  (4-3).  The  program  also  contains  a  rule  for  selecting  R  , 


the  initial  estimate  for  the  SEP.  This  rule  and  a  method  for  upgrading 


the  estimate  are  discussed  in  the  following  paragraphs, 


4.4  A  Rule  for  Selecting  R( 


As  noted  in  Sec.  4.  1,  the  computer  program  mechanizes  an  iter¬ 


ative  search  for  the  SEP.  The  search  is  started,  for  a  given  set  of 
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standard  deviations,  by  estimating  the  SEP.  If  the  first  estimate,  R0, 
causes  the  probability,  P,  to  equal  0.  5,  the  iterative  search  is  stopped 
and  the  SEP  is  equated  to  R0.  In  the  usual  case,  however,  the  first 
estimate  for  the  SEP  is  wrong  and  successive  estimates  must  be  made 
to  drive  P  toward  0.  5. 

The  first  estimate  for  the  SEP  is  made  by  setting  RQ  equal  to  the 
minimum  value  of  the  three  ratios,  °x/ 0.2,  °y/0.  2  and  °z/0.2.  This 
estimate  guarantees  that  none  of  the  normalized  standard  deviations, 
0x/Ro»  °y/ Rq  and  0z/Ro.  will  be  less  than  0.  2  during  the  first  itera¬ 
tion.  This  choice  of  RQ  avoids  the  convergence  problems  discussed 
in  Sec.  4.  2. 

If  RQ,  the  first  estimate  for  the  SEP,  is  correct,  the  computer 
will  make  one  iteration,  stop,  and  print  out  RQ  as  the  SEP.  If  the 
true  SEP,  however,  is  much  greater  than  R0,  the  computer  method 
will  not  converge  to  a  solution  for  the  SEP.  The  computer  will  make 

two  iterations,  test  for  a  diverging  solution  and  stop.  When  the  true 

I 

SEP  is  less  than  RQ,  the  computer  will  proceed  to  a  solution  for  the 
SEP  by  making  successive  estimates  for  R. 

If  the  first  estimate  for  the  SEP  is  not  correct,  the  second  esti¬ 
mate,  Rj,  is  made  by  adding  AR0,  a  small  number,  to  R0.  ARQ,  in 
turn,  is  computed  by  dividing  APQ  =  (0.5-PQ)  by  the  first  partial 
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derivative  of  P  with  respect  to  R  evaluated  at  R  =  RQ.  APQ  is  the 
error  in  probability  resulting  from  the  first  and  incorrect  estimate 
for  the  SEP.  See  Fig.  4.  1. 


Figure  4.  1  A  graph  of  P  versus  R  showing  the 

quantities  used  in  estimating  the  SEP 

The  partial  derivative,  (dP/dR)0,  gives  the  slope  of  the  probability 
curve  at  R  =  ZiQ,  A  straight  line  representation  of  this  slope  appears 
on  Fig.  4.  1.  APQ  divided  by  the  (dP/dR)Q  gives  the  approximate 
change  in  RQ  required  to  change  PQ  by  APQ. 


*  PQ  is  the  probability  calculated  for  R  =  RQ. 
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The  third  estimate  for  the  SEP  is  made  by  evaluating  thedP/dR, 
P,  AP  and  AR  at  R  =  Rj.  R^  is  computed  from  the  expression 
R2  =  Rj  +  ARj  (4-4) 

The  fourth,  fifth  and  successive  estimates  are  made  in  an  analo¬ 
gous  manner  until  the  i^  estimate  causes  the  probability  to  be  0.  5. 

The  probability,  P,  and  the  partial  derivative,  dP/£R,  are  calculated 
from  Lilliefors'  series  solution  of  Eq.  (1-1).  The  mathematical 
details  appear  in  Appendix  B. 

4.  5  Summary  of  the  Computer  Method 

The  computer  method  produces  an  exact  SEP  for  any  given  set  of 
standard  deviations.  If  the  accuracy  requirement,  however,  is  not 
too  stringent,  the  SEP  can  be  also  found  by  use  of  normalized  plots. 
These  plots,  generated  from  computer  data,  display  the  ratio  SEP/az 
as  a  function  of  the  ratio  0X/°Z.  We  have  produced  a  set  of  these  plots 
in  Sec.  5.  They  constitute  our  fourth  and  last  method  for  finding  the 
SEPi 
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SECTION  5 


THE  FOURTH  SOLUTION: 

A  GRAPHICAL  METHOD  OF  FINDING  THE  SEP 

5.  1  The  Graph 

One  SEP  can  be  computed  by  the  method  of  Sec.  4  for  any  set  of 
standard  deviations,  0  ,  and  0  .  If  a  large  number  of  SEPs  are 

X  y  Z 

calculated,  they  can  be  tabled  as  in  Fig.  5.  1. 

0X  (ft)  ov(ft)  oz(ft)  SEP  (ft) 


1 

1 

1 

1. 538 

1 

.  5 

.55 

1.045 

1 

1 

.25 

.908 

Figure  5.  1  A  sample  tabulation  of  the  SEP 
as  a  function  of  and  °z. 

The- tabulation  is  useful  for  identifying  the  SEP  associated  with  a 

i 

given  discrete  set  of  standard  deviations.  The  tabulation,  however, 

reveals  no  obvious  method  of  interpolating  for  the  SEP  associated 

with  a  set  of  values,  0  ,  0  and  that  is  not  tabled.  Much  more 

x  y  z 

information  can  be  derived  by  plotting  the  SEP  as  a  function  of  the 
standard  deviations. 
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In  Fig.  5.2  we  have  plotted  the  ratio  SEP/ versus  the  ratio 

* 

o  /o  with  o  as  a  parameter.  The  curves  are  seen  to  be  smooth 
x  2  y 

and  well-behaved.  Given  values  of  °x,  °y  and  °z>  these  curves  can 
be  used  in  a  direct  solution  for  the  SEP.  No  iterative  calculations 
are  required.  The  user  need  only  assign  values  to  the  standard  devia¬ 
tions,  form  the  ratios  0  /°  and  0  fa  ,  locate  the  appropriate  curve 

x  z  y  x 

on  the  graphs  and  read  off  the  ratio  SEP/o  .  Since  o  is  known,  the 

z  z 

SEP  can  be  computed  directly  from  the  ratio  SEP/^Z  by  multiplica¬ 
tion;  i.  e. , 

SEP  -  (-^p-)  °z  <S-» 

In  the  next  section  we  work  a  sample  problem  to  show  how  these 
curves  could  have  been  used  in  the  satellite  rendezvous  problem  of 
Sec.  1. 

5.  2  A  Sample  Problem 

I 

>In  Sec.  1  we  described  an  orbiting  laboratory  travelling  in  a 

i 

100  nm  circular  orbit.  We  said  that  the  laboratory  was  to  be  re¬ 
supplied  by  a  second  vehicle,  the  transporter,  which  would  rendez- 
vous  and  dock  with  the  laboratory.  The  rendezvous  required  a 

*  The  curves  are  patterned  on  the  CEP/oj  curves  used  by 
R.  A.  Moore  in  Ref.  3.  The  curves  were  drawn  from 
computer  data  displayed  in  Appendix  C. 


precise  prediction  of  the  laboratory's  position  at  2117  Greenwich 
Mean  Time  (Zulu).  The  prediction  included  a  SEP  of  8000  feet. 

If  tha  mission  planners  had  computed  the  SEP  by  using  the  curves 
in  Fig.  5.2,  they  would  first  have  used  three  scandard  deviations  to 
describe  the  uncertainties  in  the  laboratory's  predicted  position.  If 
the  deviations  were 


(Downrange) 

=  6350  ft. 

(5-2) 

(Crossrange) 

=  6350  ft. 

(5-3) 

(Altitude) 

=  2540  ft. 

(5-4) 

at  2117  Zulu,  the  planners,  next,  would  have  formed  the  ratios  of 
and  0  /ox. 

y 

oJa9  =  2.5  (5-6) 

x  z 

°y/°x  =  1  (5-7) 

Entering  Fig.  5.2  at  o  /a  =  2.  5,  they  would  have  proceeded  verti- 

X  z 

cally  upward  to  the  curve  for  °y/°x  =  1.  On  this  curve  they  would  have 
read,  for  the  ratio  SEP/oz,  the  value  3.  15.  Multiplying  this  value  by 
°z  =•  2540  ft.  would  have  produced  the  SEP. 

SEP  =•  3.15  x  2540  ft.  =  8000  ft.  (5-8) 

The  same  procedure  can  be  used  in  computing  any  SEP.  The 
curves  are  not  restricted  to  problems  in  orbital  mechanics.  A  SEP 


can  be  calculated  in  any  three-dimensional  uncertainty  problem  where 

2  2  2 

we  can  give  values  to  the  normal  variances  ox  *  °y  and  °z  . 
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5.  3  Application  of  the  Graphical  Solution 


In  the  graphical  method  of  finding  a  SEP,  we  assign  values  to  the 

standard  deviations  0  ,  0  and  0  .  These  values  are  generated  by  our 

x  y  z  **’ 

study  of  the  physics  and  statistics  in  the  actual  problem.  In  general, 

these  values  will  be  associated  with  a  problem  coordinate  system, 

like  Xj,  Xp,  and  x^,  and  not  with  the  coordinates  x,  y  and  z.  The 

question  arises  as  to  how  the  values  computed  for  the  three  variances 

2  2  2 

in  Xj,  x^  and  x^  space  should  be  assigned  to  ox  ,  and  oz  . 

The  solution  for  the  SEP  is  independent  of  the  order  of  values 
assigned  to  °x,  o  and  because  Eq.  (1-1)  is  symmetric.  The 
equation  is  symmetric  in  the  sense  that  the  six  combinations  of  any 
three  values  all  have  the  same  SEP.  See  Fig.  5.  3  for  one  example. 


0  x  <£t> 

0  y  ««) 

%  («) 

SEP  (ft) 

1 

2 

3 

3.  105 

1 

3 

2 

3.  105 

2 

3 

1 

3.  105 

2 

1 

3 

3.  105 

3 

1 

2 

3.  105 

3 

2 

1 

3.  105 

Figure  5.  3  A  table  showing  the  symmetry  in  the 
SEP  calculation. 


The  user's  computed  standard  deviations,  aj,  and  ay  can 

be  assigned  in  any  order  to  n  .  a  and  a  .  The  same  SEP  will 

#  x  y  z 

result  for  every  combination.  The  graphical  solution,  however,  has 

widest  application  when  the  largest  computed  deviation  is  assigned 

to  If  the  largest  standard  deviation  is  not  assigned  to  o  but  to 

0  instead,  the  ratio  0  /  0  may  be  larger  than  3.  5  and  beyond  the 

X  X  z 

curves  of  Fig.  5.2. 

Finally,  we  note  that  the  SEP  for  max,  mOy  and  is  just  m 
times  the  SEP  computed  for  0  ,  o  and  o  This  result  is  not  imme> 

a  y 

diately  obvious  from  Eq.  (1-1)  but  is  apparent  in  Fig.  5.2.  If,  for 

example,  the  standard  deviations,  a  =  5  ft,  a  =  10  ft,  and  a  = 

x  y  “ 

2.  5  ft,  are  all  doubled,  the  SEP  goes  from  3.  67  ft.  to  7.*  34  ft. 
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SECTION  6 


SUMMARY 

In  the  opening  paragraph  of  Section  1  we  advertised  four  methods 
of  solving  Eq.  (1-1)  for  the  Spherical  Error  Probable.  We  have 
described  those  four  methods  in  Sections  2  through  5.  We  have 
showed  how  a  SEP  can  be  calculated  by  approximation,  paper  and 
pencil,  computer,  or  graph. 

Our  four  methods  of  computing  the  SEP  are  not  distinctly  differ¬ 
ent.  Lilliefors'  solution  of  the  probability  integral,  for  example, 
underlies  each  of  the  methods  presented  in  Sections  3,  4  and  5.  The 
same  iterative  procedure  is  used  in  the  search  for  the  SEP.  We  have 
identified  the  methods  separately,  however,  because  different  com¬ 
puting  took  -  paper  and  pencil,  computer  or  graph  -  are  used  in 
finding  the  SEP.  The  associated  computation  times  also  differ 
markedly. 

This  paper  is  not  an  exhaustive  survey  of  methods  for  computing 
the  SEP.  There  probably  exist  more  efficient  computer  solutions  to 
the  probability  integral,  Eq.  (1  - 1 ). for  example.  We  have  only 
described,  herein,  four  methods  which  the  authors  have  found  useful 
in  making  many  repetitive  SEP  calculations  for  multiple  or  time 
varying  sets  of  standard  deviations. 


APPENDIX  A 


ESTIMATING  THE  SEP  WHEN  0  =  mo 

y  x 

A.  1  Outline  of  the  Estimation  Method 

In  Sec.  2  we  showed  how  the  SEP  (or  SEP/  °_)  could  be  estimated 

z 

in  the  special  case  where  two  of  the  standard  deviations  were  equal. 

In  this  appendix  we  develop  a  method  for  estimating  the  SEP  when 

0  =  mo  ,  m  =£  0,  1.  The  method  uses  the  SEP/o  curve  presented 

y  x  z 

in  Fig.  2.  1. 

For  purpose  of  example  we  will  only  do  calculations  for  the  case 

when  o  =  2o  .  Similar  calculations,  however,  can  be  applied  when 
y  x 

m  =  0.  5,  3,  4,  5  .  .  .  .  In  the  following  paragraphs  we  approximate 

the  SEP/  0  curve  for  0  y  *  by  calculating  the  SEP  when 
z  * 

(1)  °x  is  small  compared  to<*z. 

(2)  o  is  large  compared  to  o  . 

X  Z 

(3)  ox  equals  oz. 

i 

(4  )  <*x  equals  0.  5 <jz. 

A.  2  The  SEP  for  o„  =  2  oY  and  a  Small 

In  Sec.  2.2  we  noted  that  the  three-dimensional  probability  dis¬ 
tribution  looked  like  a  one-dimensional  distribution  when  o  and  a 

Y  x 

were  equal  and  much  smaller  thanOz>  The  same  situation  occurs 
when  o  =  2ox  and  0 x  is  small.  In  the  limit,  as  °x/o  Z'  approaches 


I 

1 


zero,  the  probability  density  function  becomes  one-dimensional  in  the 

variate  z.  The  SEP  is  given  by  the  equation: 

SEP  =0.674oz  (A-l) 

when  o  /o  =0.  In  normalized  form  the  ratio,  SEP/  o  ,  is  just  0.674. 
x  z  “ 

This  is  the  same  value  calculated  in  Sec.  2.2.  SEP/0  =  0.674  is 

z 

plotted  as  a  point  at  o  /  a  -  0  on  Fig.  A.  2. 

x  z 

A.  3  The  SEP  for  a  =  2o  and  q  Large 

y  x _ ^ _ 

For  o  =  2o  andcY  large  compared  too  ,  the  density  function, 
f(x,  y,  z)  looks  like  a  two-dimensional  distribution  in  x  and  y.  The 
SEP  is  equal  to  the  CEP  defined  in  Sec.  2.  3.  The  SEP  is  calculated 
from  Eq.  (A.  2): 


SEP  =  CEP  =  0.  59  (a  +  a  ) 

x  y 


(A -2) 


If  we  substitute  2ox  for  Oy,  add,  and  then  divide  by  oz,  the  familiar 
normalized  function,  SEP/°Z,  results: 


SEP  „  1.77  I* 
0  0 


(A  -3) 


z  z 

Equation  (A-3)  says  that  when  q  =  2o  and  ov/  o_  is  much  greater 

y  «rW  A  mi 

than  1,  the  SEP/o  z  curve  approaches  a  straight  line.  The  line  has  a 
slope  of  1.77  and  a  zero  intercept  atax/oz  =  0.  This  line  is  plotted 
on  Fig.  A.  2. 
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When  ay  =  2ox  and  ox  equals  a  z,  the  SEP  can  be  estimated  by 
using  Fig.  2.  1  and  interchanging  the  subscripts  on  the  standard  devia 
tions.  Figure  2.  1  was  plotted  under  the  assumption  that  o  =  c?x,  in 
this  new  case  we  have,  instead,  a  equal  to  o  If  y  and  z  are  inter- 
changed  wherever  they  appear  on  Fig.  2.  1,  we  get  Fig.  A.  1  and  a 
direct  solution  for  the  SEP. 


Figure  A.  1.  The  predicted  SEP/dy  curve  when  o  x  =o 
At  the  point  where  ox/Oy  =  1/2  we  read  that: 


SEP /a.. 


1.  5382 


fA-41 


Multiplying  Eq.  (A -4)  by  the  ratio  Oy  /a  z  =  2  produces  our  familiar 
normalized  function  >  SEP/o^,  and  we  have 

SEP/oz  =  2  X  1.  5382  =  2.  0764  (A-5) 

This  point  is  plotted  on  Fig.  A.  2  at  ov/a„  =  1. 

A  z> 

A.  5  The  SEP  when  o  .  =  o  =  2o 

y  z  * 

A  final  point  can  be  found  on  the  SEP  curve  for  o  =  2qv  when 

y  a 

°x  =  0.  5qz.  As  in  Sec.  A.  4  we  find  this  point  by  re-plotting  Fig.  2.  1. 
The  plot  is  not  reproduced  here.  We  interchange  the  x  and  z  sub¬ 


scripts  on  the  standard  deviations  and  find  that  the  SEP/cx  =  2.  58  at 

V°x  = 2- 

The  function,  SEP/o  x,  is  changed  to  our  standard  form  through 

multiplication  by  the  ratio  c  /o  =  1/2.  Eq.  (A -6)  results. 

X  z 


The  point,  SEP/oz  =  1.29,  is  plotted  on  Fig.  A.  2  at.  °x/°z  -  1/2. 

We  now  have  sufficient  data  to  draw  in  the  rest  of  the  SEP/o  curve. 

z 

A.  6  The  SEP/o  Curve  for  q„  =  2q 

se  y  x 


The  points  and  the  line  found  in  Sec.  A.  2  through  A.  5  appear 


xnr.n  ‘  * 


Figure  A.  2.  The  predicted  SEP/o  r  curve  for0  =  2  q 

z  y  x 

A  dotted  line  is  drawn  to  connect  the  points  and  complete  the  curve. 
The  SEP  can  be  estimated  from  this  curve  in  the  special  case  when 

O  y  -  20 

By  analogous  methods  the  function  SEP/o  z  can  be  graphed  for  any 

combination  of  o  and  o  •  i.  e.  ,  for  ov  =  mo  ;  m  =  0.  5,  2,  3,  .  .  .  . 

y  *  /  x 

When  m  =  4,  for  example,  the  approximate  SEP/o^  curve  is  generated 
by  using  the  curves  for  m  =  0,  1,2,  and  3. 


APPENDIX  B 


THE  MATHEMATICS  IN  THE  ITERATIVE 
SOLUTION  FOR  THE  SEP 


B.  1  Lilliefors1  Solution  for  P 

Our  iterative  solution  for  the  SEP  is  based  on  a  series  solution 
to  the  integral  equation: 


R  Vr2-z^ 


The  series  solution  was  reported  by  H.  W.  Lilliefors  in  Ref.  2. 

For  R  =  1  and  o  =£  o  0_,  Lilliefors'  solution  for  P  has  the  form: 

x  y  z 


*  For  °x  =  °y  =  °z  =  °»  the  probability  integral  can  be  solved  by 
a  series  expansion  of  the  integrand  in  Eq.  (A -  1 ) .  For  this  case, 
J.  S.  Toma  has  reported  that  the  SEP  equals  1.  5382  .  ^  The 
same  answer  is  derived  from  Lilliefors'  Eq.  (B-2). 
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This  solution  is  correct  for  ^  0.  If  0  =  0,  P  is  computed  from 

z  z 

the  expression: 


<-l)(n+1,<»1  <C2,n) 


(B-3) 


2n  o, 


Jn-1> 


o  n ! 
x 


The  factor,  (C2,n),  appears  in  both  Equations  (B-2)  and  (B-3). 
oj  (C2.n)  is  calculated  from  the  following  recursion  formulas: 

(C2,  n)  =  a  <*2  (C2,n)  +  b02  <C2,n) 
am  (Cz,  n)  -  1  when  m  =  n 

am  (C2,  n)  =  a«m  +  1  (C2,  n)  +  bBm+1(C2,  n)  when  m  ^  n 
Bm  (C2,  n)  =  0  when  m  =  n 

Bm  (C2.n)  =  [b«m+1(C2,n)+a6m+1(C2,n)J  when  m  =*  n 

(C2.  n)  =  0 


The  constants  a  and  b  in  the  recursion  formulas  are  defined  by 
the  equations: 

a  =  1  +  j  (k2  -1)  * 

I  b=i(k2-l) 

k  - 

k  '  °x 


*  This  expression  for  "a"  is  a  correction  to  the  expression  in 

1  2 

Ref.  2.  The  equation  in  Ref.  2  is:  a  f  1  +  1  2  (1-z  ). 

4  ov 


B.2  The  Partial  Derivative 


The  iterative  solution  for  the  SEP  requires  successive  estimates 
for  R,  the  spherical  radius.  The  computer  program  calculates  the 
successive  estimates  by  using  the  first  partial  derivative  of  P  with 
respect  to  R.  In  this  section  we  show  how  the  partial  derivative  is 
calculated  from  Lilliefors1  series  solution  for  P. 

Lilliefors'  solution,  Eq.  (B-2),  can  be  rewritten  in  the  form: 


Q  <n) 


-TSTTTT 


°x°z 


.£  N  fr-j* 

J  (°  >J 


(B  -4) 


The  standard  deviations,  0  ,  o  and  o  are  normalized  in  this 

x  y  z 

expression.  If  0  ,  0  and  o  are  redefined  to  be  the  non-normalized 

X  7  Z 


deviations,  Eq.  (B-5)  results: 


J 


/ 


p  r  ci  L 


Q  (n) 


n 


(2n-l) 


N  (n,  j) 

T 


m 


(B-5) 


We  next  create  Eq.  (B-6)  by  rearrangement  of  Eq.  (B-5). 


P  ^  C,  xp  r(2h+  q  RZj  n  (n,  j)  (B  -6) 

^  o  (2n‘1)ox  az  ^  ^ 

n  oy  x  z  j  z 

Lilliefors'  solution  is  now  in  a  form  where  we  can  take  the  deriv 


ative.  Before  taking  the  derivative,  however,  it  is  instructive  to  do 


a  partial  expansion  of  Eq.  (B-6).  The  first  three  terms  in  both  series 
appear  in  Eq.  (B-7). 


C,R3Q(1)  r  R2  R4 

p  --  0  0  -  -  N(1.0)  +  N(l.  1)  5-2  +  N<1.  2) 

y  x  z  z  z 


+  Cl^-Q(2I  N(2,0)  +  N(2,  1)5_  +  N(2,2)iL-  ... 

a  a  a  °z  i 

y  x  z  L  z  z  J 

7  2  4 

+  C}R  Q(3)  N(3,  0)  +  N(3,  1)^-5+  N(3,  2)  R  ... 

o  5  o  o  °z  o 

y  x  z  L  z  J 

•  (B-7) 

If  we  now  take  the  first  partial  derivative  of  P  with  respect  to  R,  Eq. 
(B-8)  is  formed: 


=  Cl°(‘>  3R2  N(1,0)  +  N(l.  1)S_  +  N<1.2)-* 

JR  a  °z 


0*0  a 
y  x  z 


+  ClQ(2)  SR4  N(2,  0)  +  N(2,  l)2L-r  +  N(2,2)R_ 


4-  C‘Q—  7R6  N(3,  0)  +  N(3,  1)  +  N(3,2)E 

°yb  °x  %  «* 
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+  C,Q(1)R- 
°y  °x  °Z 


2R 


N(l,l)  3N(1,2) 


+  4R 


+  C,Q(2)R 

°y3°x  °2 


N!2Jl)_+  4R3  N(2,2) 


C,Q(3)R 
°y5  °x  °z 


+  4R 

0  2 
z 


3  N(3,2) 


(B  -8) 


Equation  (B-8)  contains  four  infinite  series.  The  equation  can  be 
rewritten  into  Eq.  (B-9),  a  compact  statement  that  describes  the 
series  in  summation  form: 

dP  = 

dR 

+ 

(B-9) 


„  ^  <2n+l)  R2"  Q(n) 

1  £  % 


N(n,  j)  R 


2j 


j-0 


(2n+l) 

£  R  Q  (n) 


n-l  °y 


- • 

°x°z 


L 

J-0 


R^-'Wj. 


0, 


TT 


41 


If  we  now  multiply  and  divide  the  first  and  last  series  in  Eq.  (B-9) 
by  R,  we  produce  the  expression: 


dP  =  £1 
3R  R 


^  (2n+l)  R(2n+1)Q(n) 

(2n  -  1 ) 

n«l  o'  a  a 


oo  2  j 

N(n.j)  R 

o  ^  j 

j«0  z 


+  c 


E 

n=  1 


R<2n+1)  Q(n) 

o  " 


2j 


RZj  N(n,  j) 


(B-  10) 


Equation  (B-10)  can  be  used  to  calculate  the  <}P/dR  at  any  given 

value  of  R.  For  our  computer  application,  however,  computation 

time  is  saved  by  noting  that  Eq.  (B-10)  resembles  Eq.  (B-6).  In  fact, 

if  we  rewrite  Eq.  (B-6)  in  the  form: 

P  =  Ci  SA(n).$B(n,j)  (B-ll) 

n  T 

we  can  write  Eq.  (B-10)  as: 

=  C,  A(n)  .ZB  (n,  j) 

aR  n  R  j 

+  C,  Z  A(n)  .  Z  B(n,  j)  (B - 12) 

n  j  R 

Equations  (B-ll)  and  (B-12)  show  that  the  two  series,  ^  A(n)  and 

?B(n,  j),  are  common  to  the  solutions  for  both  P  and  the  ^P/ aR. 

J 

Given  a  o  ,  o  and  an  estimate  for  R,  both  P  and  aP/dR  can  be 
*  y  z 

calculated  from  only  one  solution  for  the  two  infinite  series. 
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APPENDIX  C 


A  COMPUTER  LISTING  WITH  DATA 


C.  1  The  SEP  Computer  Listing 


The  following  CDC  3600  FORTRAN  program  mechanizes  the  SEP 
calculations  outlined  in  Sec.  4. 


PROGRAM  SEP 

c****  this  program  computes  the  sep  for  a  probability  or  .9 

C0MM0N/S/R.SX.SY,SZ.BRiDX,DY.DZ.rR,FX,FY.F2»XN,XK.C2 
*.  I  FLAG 

COMMON/BLK/  sump  ,ALPHA(S0.60),BETA(60»80).SIGX,SJGY, 

•SI  GZ 

COMMON/X/  SIG(300)#R9(300)»TIM(300)iPR<300) 

K*  0 
PN« .  5 

C****REaD  cards 
10  CALL  READ 

N«XN 
S IZ«SZ 

IF ( EOF  >  60  ) 1000 1 12 

C******SET  RADIUS  TO  BE  M I N  POSSIBLE  VALUE  AND  YET  INSURE  CONVERGENCE 

12  R«MINlF(SX/.2.SY/.2.9Z/.2> 

I«0 

IFISY.EO.O)  GO  TO  60 

SIGX-SX/R  S  SIGYaSY/R  S  SIGZ-SZ/R  S  XK*SY/SX  S  C2*XK**2-t 
CALL  ALPM(N,C2) 

13  CALL  PR0B(N#P.SP) 

IF  (IFLAQ.NE.0)  GO  TO  14 
K*K*1 

IFtP.LE. .501.AND.P.GE. .49t.OR,X.GT.20)GO  TO  15 
R«R-(P-PN)/SP 

SIGX-SX/R  S  SIQYaSY/R  S  SIGZ-SZ/R 
GO  TO  13 

C***«  STORE  ONLY  THOSE  VALUES  WHICH  CONVERGE 

15  IMK.GT.20)  GO  TO  14 

I«  1*1 

c****»  store  values  to  be  rut  on  tare  to  be  plotted 

TIM(I)*I 
PR(I )«P 
R5( I )«R/SZ 
s  i  g  ( i >«sx/sz 

155  K«0 

sz-sz-oz 

SIGZaSZ/R 

IMSZ.LE.FZJGO  TQ  13 
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c***  write  tape 

C****  WRITE  ON  49 
16  WRITE  <10#a>#I 
WRITE(10*9) 

WRITE!  10  )  #  ( T 1  H(  K )  #  S I(l(K ) « R5(K ) » K»l#  I ) 

END  PILE  10 
PRINT  2.SX.SY 

PRINT  1* ( S IG( K ) *R5(K)iPR(K),K«l,  I ) 

2  FORMAT  < 1H1 # //30 X # *S I QX  ■  *«E20.10*  10X,*SIGY  ■  *»E20.10> 

1  FORMAT  ( 10X#  *SIQX/SIOZ*illX»  *CEP/SIQZ*. 12X,*PRO0ABILITY*i/, 

•  (3F20.8) ) 

SZ-8IZ 

SV»SY*DY 

IP(SV.LE.FY)OOTO  12 
QO  TO  10 

14  RsMINlF  (SX/.2»SY/,2iSZ/i2> 

SIQX«SX/RSSIGY«SY/R 
QO  TO  155 

60  R«MlNlF(SZ/.2»SX/,2.i. ) 

KK«0 

SIOX-SX/R 

$IGZ*SZ/R 

61  XK«SZ/SX 
C2«XK**2-1 
CALL  ALPW(N«C2) 

62  FN«1 
5Ji*-l 
PY«PPY*0 

DO  65  K«l,N 

*J1«-1*8J1 

FN«FN*K 

S!aSJl«ALPHA(l,K)/rFN*2,*«K*SlQZ**(2*K-l>*S!GX) 

S1P>II*2  *M/R 
PY«PY*SI 

IF(AB8F(SI/PY) .LT.iE-0)  QO  TO  66 

65  FPY«PPY*S!P 

66  I* (PY.QT, . 501. OR, PY.LT. ,499. AND. KK.LT.20I  GO  TO  70 
IF(KK.GE.20)  GO  TO  69 

1*1*1 

R5( I )»R/SZ 
S I G ( I ) *SX/SZ 
T I H  C I )•  I 
PR( I )«PY 
69  SZ«SZ+DZ 

IF(SZ.GT.FZ)  16*60 
70  RiR-«PY-,5)/PPY 

KK»KK*1 
S I  GX»SX/R 
SIGZ*SZ/R 
QO  TO  62 

1000  return 

6  FORMAT  (•ID**52X»*0 41 000003*13* 6X) 

9  FORMAT(*RAD!USSIG  •) 

END 
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FUNCTION  SUMJ<  I  ) 

COMMON/BLK/  SUMP  ,ALPHA(60.80),BETA(60,60>'S!QX,S!GY. 

•SIGZ 

C0MM0N/S/R(15) 

•# iflag 
Sum j*sump*o 

SJ1«-1SFJ*1 
DO  20  JJ>1«100 
JaJJ-1 
DJ»J*.5 
DO  10  K«l»! 

10  DJ* ( , 5 ) *DJ 

FJbFJ*J 

IFIFJ.EO.O. )fJ*l 
SJ1«-1*SJ1 

S  *i/(FJ*(2*S!GZ*p?)**J*DJ)*SJ1 

I F ( SUM J  .NE.O.AND.ABSFtS/SUMj  J.LT.lE-8)  GO  TO  21 

PS«S*2*J/R 

SUMP  aSUMP  *PS 

20  SUMJ  sSUMJ  *S 
1  FLAG  *i 

21  CONTINUE 

1  FORMAT  <//,<5E24.10>) 

RETURN 

END 


SUBROUTINE  PROB(N«P. IP  ) 

COMMON/BlK/  SUMP  *  ALPHA ( 80 » 80 ) » BETA< 60 . BO ) * S IGX, S I  GY* 

•  SIGZ 

C0MM0N/S/RI15) 

•« IFLAG 

DATA( TWO  PI  a. 7978845808  > 

IPaSNaO 
IFLAQaO 
SJia-1 
DO  10  I ■! # N 
SJla-1. *SJ1 

S» ALPHA ‘ 1» | )/<2.**( I*1)*SIGY«*C2*I-1)*SI0X*SIGZ)*SJ1 

ST»S*SUMJC|J 

IF ( IFLAG. NE 1 0 )  GO  TO  9 

IF ( SN . NE  •  0  •  AND . A0SF (ST/SN).LT.1E*8)G0  TOll 
SNaSN*ST 

10  SPaSP^SUMP  •S*S*( 2*I*1)/R  *ST/S 

9  IFLAQal 

11  SPaTWOP I *SP 
PaTWO  P!*SN 
RETURN 

END 
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$  ALPHA (0O»8O)»0ETA(8O*8O)#S!QX#816V» 


SUBROUTINE  ALPH(N#C2> 

COMHON/BLKs  SUMP 

*9102 
A«C2*.5^1 
B«C2*. 5 
M«N 

00  9  I  •  1  #  N 
ALPHA 

•ETA(1.I)*0 
DO  10  K*liN 
H*M-1 

DO  10  I  ■ 1  *  H 

ALPHA ( !» I*K)«A*ALPHA< 1*1. I*K)*B*BETA( 1*1* I*K) 
BETAU#I*K)»<I-1)/I*CB*ALPHAU*1.I*K)*A*BETA(  1*1#  I*K) ) 
RETURN 
END 


SUBROUTINE  READ 
C0HH0N/S/R(13) #XK#C2 
READ  1» (R( I)« I»l»13> 
PORMAT(SEIS.O) 

END 


C.  2  The  SEP  Computer  Data 


The  following  computer  data  were  used  in  plotting  the  normal¬ 
ized  curves  in  Fig.  5.2. 


SIGX  =  1.0000000000+000  SIGY  =  0.0000000000+000 


SIGX/SIGZ 

SEP/SIGZ 

PROBABILITY 

5. 00000000 

3. 52434371 

0.49933981 

3. 33333333 

2.49730070 

0.49946628 

2. 50000000 

2.01914995 

0.49985782 

2. 00000000 

1. 74079638 

0.49998594 

1 . 66666667 

1. 55608473 

0.49999961 

1.25000000 

1. 32148778 

0.49967359 

1. 00000000 

1.  17740951 

0.49999970 

0. 71428571 

1. 00505711 

0.49999860 

0. 50000000 

0.87033511 

0.49993984 

0.40000000 

0.80761956 

0.49982800 

0. 30303030 

0.  75115376 

0.49971807 

0.20000000 

0.  70566641 

0.49987442 

0. 12500000 

0.  68627760 

0.49999826 

SIGX  =  1.0000000000+000  SIGY  = 

5.0000000000-001 

SIGX/SIGZ 

SEP/SIGZ 

PROBABILITY 

5.  00000000 

4.47047233 

0.49999214 

3.  33333333 

3.08410623 

0.49971674 

2. 50000000 

2.42161704 

0.49951873 

2. 00000000 

2. 03818181 

0.49948691 

1.42857143 

1.60964950 

0.49960197 

1. 00000000 

1.28884041 

0.49970528 

0.  50000000 

0. 90782976 

0.49980035 

0. 30303030 

0.  76708754 

0.49987596 

0.  21276596 

0. 70835663 

0.49943213 

SIGX  =  1.0000000000+000  SIGY  * 

l.OOOOOOOOOOfOOO 

SIGX/SIGZ 

SEP/SIGZ 

PROBABILITY 

3.  33333333 

4.05597236 

0.49999894 

2.  50000000 

3.  12140967 

0.49988051 

2.  00000000 

2. 57805305 

0.49982290 

1 . 66666667 

2. 22537530 

0.49980750 

1.  11111111 

1.65177576 

0.49985646 

0.  50000000 

1.01956866 

0.49986146 
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0.30303030 

0.81937779 

0.49987971 

0.20408163 

0. 73849930 

0.49993855 

0. 12500000 

0.69841416 

0.49995652 

SIGX  =  1.0000000000+000  SIGY  = 

2.0000000000+000 

SIGX/SIGZ 

SEP/SIGZ 

probability 

3.33333333 

5.88984063 

0.49993903 

2.50000000 

4.47031812 

0.49996898 

2.00000000 

3.63206608 

0.49994283 

0.42857143 

2.70305127 

0.49987570 

1.00000000 

2.03910945 

0.49985054 

0.  50000000 

1.28917711 

0.49991805 

0.25000000 

0.90802420 

0.49994871 

0.20000000 

0.83356676 

0.49995688 

0.  12500000 

0.73662162 

0.49996220 

SIGX  -  1.0000000000+000  SIGY  = 

3.0000000000+000 

SIGX/SIGZ 

SEP/SIGZ 

PROBABILITY 

2.50000000 

5.84261733 

0.49921362 

2.00000000 

4. 72149463 

0.49997504 

1.66666667 

3.97779948 

0.49996857 

1.00000000 

2. 54540426 

0.49989433 

0. 50000000 

1.55325508 

0.49992435 

0.30303030 

1. 16571314 

0.49996313 

0.20000000 

0.95555704 

0.49997517 

0.14925373 

0.85103799 

0.49999301 

0.  12500000 

0.80578883 

0.49981551 

SIGX  =  1.0000000000+000  SIGY  = 

4. 0000000000+000 

SIGX/SIGZ 

SEP/SIGZ 

PROBABILITY 

2.00000000 

5.88878514 

0.49977376 

1.42857143 

4.27103567 

0.49998488 

1.00000000 

3.09053778 

0.49995657 

0. 50000000 

1.81601062 

0.49993425 

0.25000000 

1.20446927 

0.49997776 

0.20000000 

1.07690612 

0.49999859 

0. 12500000 

0.85402313 

0.50078231 

SIGX  =  1.0000000000+000  SIGY  = 

5.0000000000+000 

SIGX/SIGZ 

SEP/SIGZ 

PROBABILITY 

2.00000000 

7.  14868648 

0.49998480 
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1.42857143 

5.  14340312 

0.49999067 

1.00000000 

3.  67950015 

0.49998552 

0. 50000000 

2.  08389447 

0.49995018 

0. 30303030 

1.  50347918 

0.49997427 

0.20000000 

1. 19483273 

0. 50022996 

0. 14925373 

1. 02956778 

0.49939611 

:  =  i.oooooooooo+ooo  sigy  = 

6. 0000000000+000 

SIGX/S1GZ 

SEP/SIGZ 

PROBABILITY 

1.66666667 

7.  04731527 

0. 50021207 

1.25000000 

5.  31645805 

0.49999518 

1. 00000000 

4.  29908555 

0.49999478 

0. 50000000 

2. 36071598 

0.49996830 

0.29411765 

1. 63776678 

0.49997671 

0.  18867925 

1. 25671490 

0. 50046639 

* 
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